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Abstract 



Exciton optical absorption in self-similar aperiodic one-dimensional sys- 



^ . tems is considered, focusing our attention on Thue-Morse and Fibonacci lat- 
OO _ 

. tices as canonical examples. The absorption line shape is evaluated by solving 

t — 

f— ^ . the microscopic equations of motion of the Frenkel-exciton problem on the lat- 

On ■ tice, in which on-site energies take on two values, according to the Thue-Morse 



or Fibonacci sequences. Results are compared to those obtained in random 



1-^ ■ lattices with the same stechiometry and size. We find that aperiodic order 

a: 

O ■ causes the occurrence of well-defined characteristic features in the absorption 



spectra which clearly differ from the case of random systems, indicating a 
most peculiar exciton dynamics. We successfully explain the obtained spec- 
tra in terms of the two-center problem. This allows us to establish the origin 
of all the absorption lines by considering the self-similar aperiodic lattices as 
composed of two-center blocks, within the same spirit of the renormalization 

group ideas. 

PACS numbers: 71.35+z; 36.20.Kd; 78.90.+t 



Typeset using REVTgX 



I. INTRODUCTION 

During the last few years the notion of aperiodic order has progressively emerged to 
gain a proper understanding of new physical systems. Since the remarkable discovery of 
the quasicrystalline phasea and the technical advances in submicron physics for the fabri- 
cation of semiconductor superlattices arranged according to the (quasiperiodic) FibonacciB 
and Thue-MorseBa sequences, much work has been devoted to the study of systems whose 
structural order is described by means of deterministic substitution sequences^ leading to 
self-similar aperiodic lattices. The interest in exploring the physical properties of elementary 
excitations in one-dimensional (ID) aperiodic systems, including Fibonacci, Thue-Morse or 
Rudin-Shapiro lattices and their generalizations, goes beyond a formal theoretical analysis of 
systems deserving a simpler mathematical treatment than three-dimensional ones. In fact, 
it is actually well known that aperiodic order gives rise to novel properties which are com- 
pletely absent in both periodic (crystalline) and random (amorphous) ID systems. In this 
way, aperiodic systems exhibit highly fragmented electronETQ and phononB'El spectra that 
are Cantor sets determining the existence of critical states. These exotic electronic spectra 
strongly influence electron propagation, being somewhat intermediate between ballistic and 
difussive, which gives rise to unusual behavior of the dc conductance at finite temperature.liil 
Although these striking results were initially obtained for tight-binding and Kronig-Penney 
model Hamiltonians, we have recently shown that peculiar electronic transport properties 
must be also expected in more realistic systems.li3 

Moreover, it has also been realized that systems ordered according to the Fibonacci 
sequence exhibit some characteristic properties which are not shared by other self-similar 
aperiodic arrangements. In particular, from studies concerning the electronic spectrum 
structurejl^l Landauer resistance,LJ'Ej and phonon spectrum properties,!!^ some authors have 
claimed that the kind of order associated to the Thue-Morse sequence must be considered 
as intermediate between the quasiperiodic order displayed by Fibonacci systems and the 
usual periodic order. We believe that this result is not surprising. In fact, the Fourier 



spectrum of Fibonacci lattices is pure point, hence indicating the existence of long-range 
order associated to the quasiperiodic nature of the underlying lattice, whereas the Fourier 
spectrum of Thue-Morse lattices is singular continuous. Ii3 Nevertheless, we feel that, in 
some of the above mentioned works, the criteria introduced to determine the degree of 
periodicity associated to a given self-similar lattice are somewhat vague and could lead to 
possible misinterpretations.Ej Therefore, to attain a deeper insight into the nature of the 
order displayed by different kinds of aperiodic structures, it seems convenient to investigate 
transport properties different from those usually considered (electron propagation, phonon 
dynamics) . 

In this regard, optical properties of aperiodic lattices have received much less attention 
and, to our knowledge, most of the work has been restricted to the study of optical phe- 
nomena in Fibonacci superlattices.cJ'Ey In this work we investigate the optical absorption 
spectra of two different kinds of self-similar aperiodic systems, namely the Fibonacci lattice 
(FL) and the Thue-Morse lattice (TML), and compare them with the optical spectra char- 
acteristic of both binary random and periodic related systems. To this end, we make use 
of a general treatment which allows us to study the dynamics of Frenkel-excitons in these 
lattices, solve the microscopic equations of motion and find the optical absorption spectra. 
This study is inspired in our previous work showing that short-range correlated disorder has 
profound effects on trappingtil and optical propertiesc3 of Frenkel-exciton systems. These 
results lead, in a natural way, to the question as to whether long-range aperiodic order mod- 
ifies exciton dynamics in comparison to long-range disorder effects. The main conclusions of 
this work are twofold. First, we show that both FLs and TMLs exhibit optical absorption 
spectra quite different from those obtained in random and periodic lattices. Therefore, op- 
tical spectra can be used to characterize experimentally the occurrence of aperiodic order in 
the sample. Second, we show that optical spectra are able to discriminate also the particu- 
lar kind of aperiodic order present in the system. Hence, the analysis of optical absorption 
spectra appears as an excellent diagnostic tool to characterize the structural order from an 
experimental point of view. 



The remaining of the paper is organized as follows. In Sec. II we describe our model and 
the different self-similar aperiodic arrangements we are going to investigate, and we show 
how optical spectra can be numerically obtained. In Sec. Ill we give a detailed account of 
the main lines appearing in the spectra, and we compare them with spectra of random and 
periodic systems. Section IV is devoted to find a relationship between those features and 
the underlying lattice topology by means of the so-called two-center problem, guided by 
renormalization group ideas. Then, the origin of the main lines appearing in the spectra is 
explained in Sec. V on the basis of the two-center model. Section VI concludes the paper with 
some general remarks on the physical implications and possible extensions of our results. 

II. MODEL 

We consider a system of A^ optically active centers, occupying positions in a regular 
ID lattice with spacing unity. For our present purposes we neglect all thermal degrees of 
freedom, and thus we omit electron-phonon coupling and local lattice distortions. Therefore, 
the effective Hamiltonian that describes the Frenkel-exciton problem can be written in the 
well-known tight-binding form with nearest-neighbor interactions as follows (we use units 
such that h = 1) 

T-(- = Y^ Vkalak + T ^ (4«fc+i + 4+i«fc)- (1) 

k k 

Here a], and a^ are Bose operators creating and annihilating an electronic excitation of 
energy Vk at site k, respectively. T is the nearest-neighbor coupling, which is assumed to 
be constant in the whole lattice. In what follows we consider that Vk can only take on two 
values, Va and Vb, and we shall arrange them either aperiodically, according to the Thue- 
Morse and Fibonacci sequences, or randomly. For convenience we define ca {cb = 1 — ca) as 
the ratio between the number of sites A (B) and the total number of sites N in the lattice. 
TML and FL are canonical examples of deterministic and aperiodically ordered systems, 
and they can be generated by the following substitution rules: A -^ AB, B -^ BA for the 



Thue-Morse sequence and A -^ AB, B ^ A for the Fibonacci one. In this way, finite and 
self-similar aperiodic lattices are obtained by n successive applications of the substitution 
rule. The nth generation lattice will have 2" elements in the TML and F„ elements for 
the FL, where F„ denotes a Fibonacci number. Such numbers are generated from the 
recurrence law F^ = i^n-i + -^n-2 starting with Fq = Fi = 1; as n increases the ratio 
Fn~i/ Fn converges toward r = {\/h — l)/2 = 0.618 . . ., an irrational number which is known 
as the inverse golden mean. Therefore, the on-site excitation energies are arranged according 
to the sequence VaVbVbVaVbVaVaVb ■ ■ ■ in the TML and VaVbVaVaVbVaVbVa- ■ ■ in 
the FL. The values of ca and cb are strictly equal to 0.5 for any generation of the TML. On 
the contrary, the values of ca and cb depend on the particular generation of the FL, but for 
large systems one has ca ^ t and c^ ~ 1 — t. Finally, it is worth noticing that i?-centers 
appear isolated in FLs. This is an important fact in order to explain the results we will 
present later. 

Having presented our model we now briefiy describe the method we have used to calculate 
the absorption spectra. The line shape I{E) of an optical-absorption process in which a single 
exciton is created in a lattice with A^ sites can be obtained as follows." Let us introduce a 
set of correlation functions 

Gfc(t) = E(0|«feW«]|0)' (2) 

i 

where |0) denotes the exciton vacuum state and ak{t) = exp(z7it)afcexp(— -jTit) is the anni- 
hilation operator in the Heisenberg representation. The function Gk{t) obeys the equation 
of motion 

tj^G,{t) = J2H,,G,{t), (3) 

with the initial condition Gfc(O) = 1. The diagonal elements of the tridiagonal matrix H^j 
are Vk whereas off-diagonal elements are simply given by T. The micrsocopic equation of 
motion is a discrete Schrodinger-like equation on a lattice and standard numerical techniques 
may be applied to obtain the solution. Once these equations of motion are solved, the line 
shape is found from the following expression 



I{E) = -^j^ dt e-"* sm{Et) Im (j^ Gk{t)] , (4) 

where the factor exp(— at) takes into account the broadening due to the Lorentzian instru- 
mental resolution function of width a. 



III. RESULTS 

We have solved numerically the equation of motion (^ using an implicit (Crank- 
Nicholson) integration scheme. In the remainder of the paper, energy will be measured 
in units of T whereas time will be expressed in units of T~^. Aperiodic lattices are gen- 
erated using the inflation rules discussed above. We have checked that the main features 
of the spectra are independent of the system size. Henceafter we will fix A^ = 2^^ = 2048 
for the TML and A^ = Fiq = 1597 for the FL as representative values. In addition, stan- 
dard random generators are used to obtain disordered lattices with the required size A^ and 
value of cb {N = 2048, cb = 0.5 to compare with TML and A^ = 1597 and cb = 0.382 to 
compare with FL). In order to minimize end effects, spatial periodic boundary conditions 
are introduced in all cases. Once the functions Gk{t) are known, the line shape I{E) is 
evaluated by means of (^. Since we are mainly interested in characterizing the effects due 
to aperiodic order as compared to randomness rather than in a detailed discussion of the 
optical absorption process, we will fix the values of Va, Vb and T, focusing our attention on 
the comparison between different types of arrangements of optical centers. Furthermore, in 
order to facilitate the comparison with our previous work, we have set V^ = 4, Vg = 10 and 
T = —1 henceafter. The width of the instrumental resolution was a = 0.5. The maximum 
integration time and the integration time step were 16 and 8 x 10""^, respectively; larger 
maximum integration times or smaller time steps led to the same general results. 

For the sake of clarity let us consider, in the first place, the typical spectra associated 
to both pure A and pure B lattices corresponding to periodic cases. In the pure A lattice 
the spectrum is a single Lorentzian line centered at -E^^g = Va + 2T, which with our choice 



of parameters is -E^^e = 2.0. When _B-centers are introduced at random in the lattice, 
a broadening of this main hne is observed accompanied by a shift of its position towards 
higher energies. In random systems both the broadening and the shift increase on increasing 
defect concentration,c3 in agreement with the average- T-matrix approximation (ATA).c3 A 
similar behavior takes place in the pure B lattice when A-centers are introduced with the 
main difference that, in this case, the single Lorentzian absorption line is originally located 
at -Ej^re = Vb + 2T = 8.0. However, ATA is no longer valid to determine optical spectra 
in TML and FL due to the long-range correlation induced by aperiodic order, as we will 
discuss later. 

Keeping these general results in mind we now proceed to discuss the main features of 
the spectra obtained in aperiodic systems. We shall start with the TML and compare it to 
a typical random lattice. The obtained results are shown in Fig. |I|, where all the spectra 
have the same area. From a close inspection of this figure several conclusions can be drawn. 
First of all, we observe the occurrence of a strong line centered at E = 2.9 in the TML. This 
line is accompanied by a small shoulder at around E = 3.8 (the position of the shoulder 
has been obtained using two Lorentzian functions to fit data in the energy range from 
up to 6). Moreover, two satellites appear in the high-energy region of the spectrum at 
about E = 9.0 and E = 10.2. On the other hand, concerning the random lattice, we 
note that the main absorption line is centered at about E = 2.6, closer to the position 
corresponding to the single line in the pure A lattice although the small shoulder remains 
at about E = 3.8. In addition, the intensity of the overall absorption features in the energy 
range < i? < 6 is smaller than those corresponding to the TML. To conclude, we observe 
that the random lattice also presents a characteristic pair of satellites at the high-energy 
region of the spectrum. One of them is centered aX E = 10.2, as occurs in the TML, but the 
other is found at E = 8.2. Finally, the satellite at ii^ = 9.0, clearly observed in the TML, 
appears as an almost unnoticeable shoulder causing the asymmetry of the line at ii^ = 8.2. 

Let us now turn to the FL and compare it to a random lattice with the same size N 
and cb- Results are shown in Fig. |^. The FL presents two clearly distinct lines. In the 



low-energy range we observe a main absorption line centered E = 2.9 embodying an almost 
unobservable shoulder ai E = 3.8 meanwhile, in the high-energy region of the spectrum, a 
single satellite at E = 10.2 is observed. On the other side, the spectrum associated to the 
random lattice shows a main absorption line at about E = 2.3 along with a smaller shoulder 
at about E = 3.8 in the low-energy region whereas, at higher energies, two broad satellites 
are clearly observable at E = 9.0 and E = 10.2. 

By comparing Figs. |l| and ^ we are led to the conclusion that, on the basis of the observed 
optical absorption spectra, very significant differences exist, not only between aperiodic sys- 
tems and the corresponding random ones, but also between the two realizations of aperiodic 
order we have considered. In fact, on the one hand, besides the small shoulder above men- 
tioned, there are only two distinct lines in the absorption spectrum of the FL whereas three 
different lines are clearly observable in the spectrum corresponding to the TML. On the 
other hand, the absorption lines at E = 2.9 and E = 10.2 are more intense in the FL 
spectrum than in the TML one. Finally, the satellite peak at ii^ = 9.0, clearly visible in the 
TML spectrum is not observed at all in the FL case. In the following Section we explain 
the origin of these characteristic features in terms of short-range quantum effects. 

IV. THE TWO-CENTER PROBLEM 

One of the most remarkable aspects of the electronic spectra in ID aperiodic systems is 
their highly fragmented nature which corresponds to a Cantor-like set with zero Lebesgue 
measure in the thermodynamical limit. The fragmentation pattern of the spectra varies on 
increasing the generation of the lattice and their detailed structure is mainly determined 
by short-range effects. This point was earlier suggested by means of the real space renor- 
malization group, where the number of energy levels appearing at the first stage of the 
renormalization process determines the number of main clusters in the spectra.l^sE^ This 
number of levels depends on the adopted blocking scheme which, for a binary system within 
the weak bound approach, usually decouples the original lattice in a series of single {A, B) 
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or double {AA, AB, BB) constituent elements. The procedure just outlined justifies the 
purported asymptotic stability of the electronic spectra of electronic Fibonacci systems.liil 
Furthermore, we have recently shown that the great success of the renormalization group 
scheme can be directly traced back to the fact that systems generated from the application 
of a substitution sequence encode more information, in the Shannon sense, than classical 
periodic systems.c3 As a consequence, from the assumption that the lattice topology must 
have profound influences on the exciton dynamics, it seems natural to extend the main 
ideas inspiring the renormalization procedure to account for the origin of the different lines 
appearing in the optical spectra of aperiodic systems. 

To this end, we shall consider the two-center problem describing the optical absorption 
spectrum of two isolated but coupled sites, labelled 1 and 2. 

ijGi{t) = ViGi{t) +TG2{t), (5a) 

i^G2{t) = V2G2{t) +TGr{t). (5b) 

at 

Solving these equations exactly with the initial conditions 6*1(0) = ^2(0) = 1 and inserting 
the result in (^ one obtains (we neglect the instrumental resolution for simplicity) 

I{E) = I+6{E-E+) + I_5{E-E^), (6) 



where 



2^2 



V l^ ^ 21 -1/2 

1 2T 



(7) 



and 



As expected, the optical absorption spectrum of the two-center problem presents two 
well-defined lines. From the Eqs. (^ and (|]) we arrive at the following possible situations. 
If Vi=V2, the intensity /+ vanishes so that spectrum exhibits a single line. Depending on 
the nature of the centers this line will be centered either at E^^ = Va+T = 3.0 {AA pairs) 

9 



or E^^ = Vb + T = 9.0 {BB pairs). On the contrary, if the on-site excitation energies are 
different (say V\ = Va and V2 = Vb) the optical spectrum presents two components centered 



at E^^ = {Va + Vb)/2tT^1 + {Va - VbY/4T^ = 7± v^ so that, in our units, E^^ ~ 3.8 
and E^^ ~ 10.2. 

Therefore, with the aid of the two-center model we can uniquely assign specific absorption 
lines to each of the pairs in which our original lattice can be decomposed, according to the 
renormalization group ideas mentioned above. In this sense, the signatures of AA and BB 
pairs are single lines located at E^^ = 3.0 and E^^ = 9.0. In addition, AB or BA pairs 
can be associated, irrespectively, to the simultaneous presence of two characteristic lines 
in the spectrum, centered at E^^ ^ 3.8 and E^^ ~ 10.2. According with these precise 
assignments, the origin of the main lines and satellites appearing in the absorption spectra 
of both FLs and TMLs can be unequivocally established. 

V. DISCUSSION 

In this section we explain the origin of the lines appearing in Figs. |l| and |^ making use 
of the two-center results. For convenience, we shall discuss both kinds of aperiodic lattices 
separately. 

A. Thue-Morse absorption spectrum 

Let us focus our attention on Fig. |I|. The main line centered at E = 2.9 is very close to 
the characteristic line E^"^ = 3.0 associated to the AA pair, hence strongly suggesting the 
possible origin of this absorption line. At this point it is important to note that the ATA 
approach cannot account for the presence of this line, since the shift of the main Lorentzian 
at -E^j-e = 2.0 in the pure A lattice due to the presence of _B-centers in a concentration 
cb = 0.5 amounts to only 0.6 units. This value is more than 30% lower than that obtained 
in the TML spectrum. Thus, it becomes clear that the main absorption line observed in the 
TML is not simply the -E^^e = 2.0 line shifted by the presence of -B-centers, as occurs in 
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the random lattice. This result suggests that the aperiodic order displayed by the TML has 
profound effects on the resulting exciton dynamics, which in turn manifest in the optical 
spectra. To find a heuristic explanation of the different exciton dynamics in TML, we would 
like to draw the attention on the fact that the A-centers can appear only isolated or grouped 
in pairs in the TML, but never forming larger groups, as it could be the case in random 
lattices. In fact, the presence of these larger segments, which behave locally as pure A lattice 
segments having _B-centers at the ends, is what explains the shift of -E^^e = 2.0 line towards 
higher energies within the framework of ATA. Therefore, the absence of large groups of 
y4-centers in the TML, along with the relative abundance of AA pairs instead, causes the 
occurrence of a noticeable and high peak at ~ 3.0. Similar reasoning explains the absence 
of the -E£jre ~ ^-0 li^^ i^ ^^^ TML spectrum, whereas such a line is clearly seen in the 
corresponding random lattice spectrum, shifted to £" = 8.2 due to the presence of A-centers. 
Furthermore, the marked satellite at E^^ = 9.0 must be associated to the presence of many 
BB pairs in the TML. Finally, the characteristic absorption satellites associated to AB pairs 
are revealed as a shoulder {E^^ = 3.8) at the high-energy side of the main line and as a 
absorption line at E^^ = 10.2. In order to further confirm this identification we have made 
use of Lorentzian fitting of data to evaluate the ratio between the relative intensity of lines 
at E^^ = 3.8 and E^^ = 10.2, which is found to be I-/I+ — 2.04. This value agrees rather 
well with the theoretical estimation I-/I+ = {VTO + l)/(vTO — 1) ~ 1.92 obtained from 
Eq. ®. 

B. Fibonacci absorption spectrum 

Now we turn our attention on the Fig. 0. Once again the contribution due to AA pairs 
(main peak at 2.9) and AB pairs (small shoulder at ~ 3.8 and satellite peak at 10.2) are 
clearly seen in the absorption spectrum, hence supporting the convenience of our the two- 
center description. Moreover, one of the most remarkable characteristic of this spectrum, 
as compared to that corresponding to the TML, is the dramatic absence of the E^^ = 9.0 
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line. In fact, according to our previous discussion such a line comes from the contribution 
of BB pairs but, as it is well known, such pairs are forbidden in the FL. We feel this is a 
very significant result since it further confirms the correctnes of our interpretation about 
the origins of the different lines appearing in the spectra and, at the same time, allows for 
an easy and confident differentiaton between different kinds of aperiodic self-similar lattices 
from an experimental point of view. 

VI. CONCLUSIONS 

In summary, we have studied the absorption spectra corresponding to the Frenkel-exciton 
Hamiltonian on self-similar aperiodic systems described by the Thue-Morse and Fibonacci 
sequences. By comparing the obtained spectra with those corresponding to random lattices 
we conclude that FLs and TMLs exhibit characteristic absorption spectra, different in many 
aspects from those of binary random lattices with the same stechiometry, and that certain 
spectral lines can be used to characterize the aperiodic order associated to FLs from that 
related to TMLs. On the other side, from the viewpoint of physical applications, we have 
obtained analytical expressions which explain our spectra and relate microscopic system 
parameters like on-site excitation energies, to experimental data like position and strengths 
of the lines. This relationship surely should facilitate future experimental work on optical 
properties of quasicrystalline solids. 

Our treatment allows us to introduce, in a rather straightforward and natural way, con- 
cepts inspired in renormalizaton group techniques, which have accomplished a great success 
in describing the electronic spectra of aperiodic systems. On the light of the obtained results 
and previous discussions, we think that the question as to whether Thue-Morse systems are 
more or less periodic than Fibonacci ones, a controversy which has raised some debate during 
the last few years, is still ill posed. In our opinion, both Thue-Morse and Fibonacci systems 
display a new kind of order, namely self-similar aperiodic order, which has its own peculiari- 
ties, and cannot be compared with periodically ordered systems in a simple way. In fact, we 
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have recently shown that self-similar aperiodically ordered systems are able to encode more 
information, in the Shannon sense, than usual periodic onesCJ, thus opening a new way to 
aswer this question. This line of reasonings may lead to a novel vision on the concept of 
order. Rather than to think into different kinds of order, classified into separate categories 
which are compared in a quantitative way (in the sense above mentioned of a particular 
category to be less random or more periodic than any other one), maybe more fruitful to 
think into different hierarchies of order. This perspective, which is inspired into the mathe- 
matical relationships between periodic, quasiperiodic and almost periodic functions, might 
be of interest to those researchers working on this field. 
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FIGURES 
FIG. 1. Absorption spectra for a ID Thue-Morse lattice (solid line) and a random lattice 

(dashed line). In both cases the system size is A^ = 2048 and the concentration of S-centers is 

CB = 0.5. 

FIG. 2. Absorption spectra for a ID Fibonacci lattice (solid line) and a random lattice (dashed 
line). In both cases the system size is A = 1597 and the concentration of i?-centers is cb = 0.382. 
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